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Представлена модель энергосистемы, работающая в реальном масштабе времени с неограниченным интервалом моделирова- 
ния, предназначенная для исследования быстрых электромагнитных и длительных электромеханических процессов. В состав 
модели входят два тепловых энергоблока с синхронными генераторами, синхронные и асинхронные машины, источник неогра- 
ниченной мощности, блочные трансформаторы, линии электропередачи, пассивная нагрузка и коммутационный блок для ими- 
тации повреждений типа короткое замыкание и разрыв. Показаны соответствие параметров модели паспортным характеристи- 
кам синхронной машины, верное поведение блока в пусковых, переходных и стационарных режимах. Описаны возможности 
модели, язык формирования заданий и средства визуализации результатов. 


В настоящее время происходит создание малых 
технологических электрических систем на базе мест- 
ных или сопутствующих первичных энергетических 
источников. В нефтегазовой отрасли распространя- 
ются тепловые газотурбинные станции. Мощности 
таких систем зачастую покрывают полностью соб- 
ственные производственные нужды и тяготеют к ав- 
тономной работе. Распространенный в больших 
электроэнергетических системах (ЭЭС) способ 
управления основан на разделении движений на 
медленные и быстрые [1,2]. Эта технология оправда- 
на и оптимальна для крупных ЭЭС. Для анализа со- 
ответствующих режимов и процессов создан специа- 
лизированный аналитический и инструментальный 
аппарат, эксплуатация которого осуществляется от- 
делами узких специалистов. Используемые при этом 
модели «подтянуты» к действующей концепции раз- 
деления единого процесса на несколько и, отражая 
главные свойства ЭЭС для соответствующего вычле- 
нения, каждая отдельная модель теряет некоторую 
часть свойств реальной системы, т.е. огрубляет дей- 
ствительное поведение объектов. Использование уз- 
коспециализированного инструментария решает ос- 
новную долю эксплуатационных задач. Однако, в 
случаях, когда развитие аварийных ситуаций охваты- 
вает широкий спектр временных возмущений - от 
быстрых электромагнитных до медленных систем- 
ных движений, возникает потребность интегрирова- 
ния разрозненных результатов в целостную картину. 
Здесь приходится уповать на опыт и искусность ис- 
следователей. Применение полнофазных физически 
подобных моделей устраняет описанный методиче- 
ский недостаток технологии анализа поведения си- 
стемы по частям [3]. Положительным свойством эт- 
их моделей является возможность непосредственно- 
го сопоставления результатов моделирования с запи- 
сями регистраторов аварийных ситуаций. Такие мо- 
дели максимально приближены к свойствам реаль- 
ных энергосистем. Программирование и восприятие 
результатов не требует существенного абстрагирова- 
ния от реальности и поэтому, кроме исследователь- 
ских функций, они могут выполнять функции трена- 
жеров работников различных специальностей, а не 
только диспетчерского персонала. 


Математическая имитационная модель с позиции 
исследователя воспроизводит определенные свой- 
ства ЭЭС точно также, как и физическая модель. Ци- 
фровые, аналоговые и физические модели мало при- 
годны для получения результатов в аналитической 
форме. В то же время они позволяют изучать реаль- 
ный объект, в нашем случае - ЭЭС, оптимизировать 
его поведение, получать качественные выводы о воз- 
можности функционирования даже в аварийных ре- 
жимах. И все это может быть многократно повторено 
в ходе экспериментов с разными условиями. 

Между физической и математической имита- 
ционной моделью есть различия [4]. Физически по- 
добные модели чаще всего работают в реальном 
времени или в ускоренном. Основное требование 
физического воплощения - необходимая степень 
подобия. Закономерности функционирования ре- 
ального объекта в них воспроизводятся автоматиче- 
ски. Отпадает необходимость выявления и матема- 
тической формализации этих закономерностей, по 
сути, отпадает необходимость познания способа 
функционирования элементов объекта. В компью- 
терных или аналоговых моделях основная нагрузка 
по обеспечению адекватности ее поведения поведе- 
нию реальному объекту падает на математически 
формализованные знания о свойствах объекта в це- 
лом и его элементов. Аналоговые электрические 
модели ЭЭС [3], подобно физическим, также обес- 
печивают высокую скорость моделирования. Их 
вычислительная эффективность практически про- 
порциональна сложности математических моделей, 
т.к. время решения не зависит от размерности си- 
стемы дифференциальных уравнений. Однако у ап- 
паратных аналоговых вычислителей есть суще- 
ственные недостатки: высокая стоимость моделей, 
отсутствие универсализма (узкая специализация), 
принципиальная невозможность сохранения со- 
стояний, что замедляет проведение экспериментов. 

Компьютерные модели свободны от перечи- 
сленных недостатков аналоговых вычислителей. В 
их основе, как и у аналоговых, математическая мо- 
дель в форме нелинейных алгебро-дифферен- 
циальных уравнений, состоящая, в свою очередь, 
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из моделей компонентов. Особые составляющие 
компьютерной имитационной модели - это чи- 
сленные схемы, соответствующие используемым 
методам интегрирования, программа и выполняю- 
щее ее аппаратное обеспечение, в простейшем слу- 
чае - персональный компьютер. 



Рис. 1. Взаимосвязь методов, математических методов и ха- 
рактеристик имитационной модели 


Основными качественными характеристиками 
имитационной модели является точность и ресур- 
соемкость, которую можно выразить через затраты 
машинного времени на шаг интегрирования. Вза- 
имосвязи, показанные на рис. 1, на практике озна- 
чают, что попытка существенно улучшить что-ли- 
бо, например, добавить дополнительные уравне- 
ния в модели, скорее всего, негативно отразится на 
остальном - в частности, на скорости выполнения 
расчетов и точности. Кроме того, модели должны 
соответствовать методам, которые в свою очередь, 
должны гарантировать устойчивость расчетов. 
Особенность стандартных пакетов состоит в том, 
что заранее неизвестны ни модели, ни окружение 
моделирования, поэтому при их разработке исхо- 
дят из худшего - стараясь в таких условиях гаран- 
тировать “надежность”, порою закрывая глаза на 
скорость моделирования и даже на точность. В то 
же время, отталкиваясь от конкретных целей моде- 
лирования, опираясь на свойства конкретной мо- 
дели, понимая, где допустимы компромиссы, что 
известно только исследователю, и, порою, выявля- 
ется только в ходе экспериментов, можно повы- 
шать качество моделирования многократно [5]. Ва- 
риант созданной имитационной модели ЭЭС 
(рис. 2), является результатом множества таких 
компромиссов и сбалансирован в отношении ско- 
рости и точности. Относительно медленные про- 
цессы в парогенераторе с постоянными времени на 
уровне десятых долей секунд и более моделируют- 
ся явными методами; относительно быстрые про- 
цессы в электрических линиях с постоянными в де- 
сятки миллисекунд и менее - неявными. Ядро мо- 
дели написано на языке С++. Функции моделиро- 
вания работают в отдельных потоках. Всё это по- 
зволило достичь скорости расчетов 50 мкс на шаг 
интегрирования при использовании процессора 
АіЫоп ХР 2500 и числе уравнений около двухсот, 
что соответствует режиму моделирования в реаль- 
ном времени. При моделировании ЭЭС на частоте 


50 Гц шаг 500 мкс уже обеспечивает достаточную 
точность в неаварийных режимах. При исследова- 
нии быстрых процессов при коротких замыканиях 
шаг может составлять единицы микросекунд и ме- 
нее. Использованные приемы программирования 
позволили применить высокоточные модели син- 
хронных машин и тепловых энергетических уста- 
новок с системами регулирования и автоматики, а 
также модели трансформаторов, линий электропе- 
редач и нагрузок [3]. Схема ЭЭС, с помощью кото- 
рой производилась идентификация свойств эл- 
ементов и моделирование процессов, приведена на 
рис. 2. 
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Рис. 2. Структурная схема цифровой модели: С - система ; 
ТГ1, ТГ2 - турбогенераторы ; Т1, Т2 - трансформаторы; 
ЛС, ЛІ, Л2 - линии электропередач, БК ~ коммута- 
ционный блок; ВС, ВГ и ВЛ ~ выключатели 


Первичная силовая установка представлена ис- 
точником рабочего тела (пара или газа) с програм- 
мно-управляемыми параметрами и турбиной с ци- 
линдрами высокого, среднего и низкого давлений с 
системой регулирования. Система автоматическо- 
го регулирования турбины (АРТ) представлена че- 
тырьмя элементами. Во-первых, автоматическим 
регулятором скорости вращения турбогенератора с 
зоной нечувствительности около 3...4 % и устрой- 
ством коррекции начальной неравномерности, 
компенсирующим влияние на приемистость тур- 
бины пара, аккумулированного в промперегрева- 
телях. Во-вторых, автоматикой предотвращения 
повышения скорости вращения турбины до уста- 
новки срабатывания центробежных выключателей 
турбины при быстром сбросе нагрузки генератора. 
В-третьих, точным регулятором мощности турби- 
ны, обеспечивающим вторичное регулирование 
частоты и активной мощности. И, наконец, кана- 
лом быстрой кратковременной разгрузки турбины 
и быстрого длительного ограничения мощности по 
сигналам противоаварийной и режимной автома- 
тики энергосистемы. Реальным прототипом моде- 
ли АРТ является электрогидравлическая АРТ, элек- 
трическая и гидравлическая части которой увязы- 
ваются в электрогидравлическом преобразователе. 
Датчики системы АРТ представлены инерционны- 
ми звеньями, по свойствам подобными механиче- 
скому и электрическому датчикам скорости, элек- 
трическим датчикам активной мощности генерато- 
ра и давления пара в линии промперегрева. Испол- 
нительными элементами АРТ являются гидравли- 
ческие сервомоторы регулирующих клапанов ци- 
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Рис. 3. Результат моделирования процесса синхронизации СГ без нагрузки при подаче номинального напряжения возбужде- 
ния на обмотке 


линдров высокого и среднего давлений и предста- 
влены в модели соответствующими инерционны- 
ми звеньями с ограничениями [7, 8]. 

Математическая модель синхронного генерато- 
ра использована в форме системы дифференциаль- 
ных уравнений в ортогональных координатах (Ц, 
жестко связанных с его ротором, с представлением 
параметров статора и ротора во взаимной системе 
относительных единиц [1—3]. В модели заданны 
следующие пространственные соотношения: оси 
(ІЧ вращаются против часовой стрелки; ось (1, сов- 
падающая по положению с осью ротора, опережает 
ось ч. В описываемой модели не учтены потери в 
стали и ее насыщение, наличие емкостей между об- 
мотками и прочими деталями внутри машины, а 
также зависимость активных сопротивлений обмо- 
ток от температуры. Ротор машины представлен 
многоконтурной схемой замещения в виде нес- 
кольких параллельно включенных активно-индук- 
тивных цепочек с постоянными параметрами в 
двух осях й и ч [3]. Переход от токов и напряжений 
модели синхронной машины (СМ) в осях <Ц к то- 
кам и напряжениям трехфазной система АВСО вы- 
полнен в соответствии с преобразованиями Парка- 
Горева [1—3]. 

Для визуализации и обработки параметров ис- 
пользованы библиотеки из пакета А8ІМЕС [6], ко- 
торый создается на кафедре промышленной элек- 
троники Томского университета систем управле- 
ния и радоэлектроники. Программа позволяет на- 
страивать параметры моделирования, наблюдать 
значения переменных и управлять выводом графи- 
ков посредством меню, консоли и скриптообра- 
зующего языка. С помощью применения техноло- 
гии многопоточности, все операции можно выпол- 
нять во время моделирования. По ходу работы ве- 
дется протокол команд. Наиболее употребитель- 
ные команды: кеі - устанавливет значения параме- 


тров модели; $іт N - запускает моделирование N 
шагов; $аѵе и Іоасі - сохраняют и загружают состоя- 
ния. Пример: 

зеС ЕЬА 30; // приравнять КЬА (активное со- 
противления нагрузки в фазе А) 30; 

зіт 1000; // выполнить 1000 расчетов с те- 
кущим значением шага интегрирования 

заѵе "Гііе"; // сохранить текущее состоя- 
ние модели в файл Еііе.зср 

ІоаФ "Гііе"; // загрузить данные из Іііе.зср 

Общее представление о свойствах и функцио- 
нировании трехфазной цифровой модели можно 
получить через описание экспериментов. На рис. 3 
приведена запись эксперимента по синхронизации 
ненагруженного генератора ТГ1 с сетью. Невоз- 
бужденная явнополюсная машина с включенной 
системой АРТ и отключенным автоматическим ре- 
гулятором возбуждения (АРВ) раскручена до ско- 
рости, близкой к синхронной скорости системы; 
ТГ подключается к системе в произвольный мо- 
мент времени (внутренний взаимный угол йеііа бо- 
лее +90°) с заметными возмущениями внутреннего 
состояния; на 12-й секунде скачком подается но- 
минальное возбуждение Ш= 1 , растут ток обмотки 
возбуждения К, внутренняя ЭДС Ещ и снижаются 
практически до нуля угол йеііа и внутренняя ЭДС 
ЕМ. Наблюдаемое поведение соответствует поведе- 
нию реальных СМ. 

На рис. 4 изображен следующий за синхрониза- 
цией ТГ этап быстрой загрузки до номинальной 
мощности включенного в систему турбогенератора 
при работающих АРТ и АРВ. С целью выявления 
скоростных возможностей теплосиловой установ- 
ки условия опыта изменены относительно суще- 
ствующих на практике, для чего постоянная време- 
ни механизма управления турбиной снижена до 
ТпшІ= 1 с. Ускоренная подача пара на турбину воз- 
мущает систему регулирования турбины. Заметны 
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значительные низкочастотные колебания, обусло- 
вленные АРТ. Однако нарушения устойчивости в 
целом не происходит. Внутренний взаимный угол 
(іеііа увеличивается и стабилизируется на уровне 
более 90°. 

Для идентификации свойств модели в области 
быстрых электромагнитных процессов проведены 
два эксперимента по имитации однофазного и 
двухфазного коротких замыканий вблизи шин 
блочного трансформатора. Параметры СМ мощно- 
стью 67 МВА - Шом=10,5 кВ; Іном=3,68 кА; 
Хаб=0,85; Хад=0,45; Х$і§та=0,15; Хы§та]3=0,18; 
Х$і§таС>=0,23; Хы§та{=1,03 - взяты из примера 
9.1 18], а параметры трансформатора 10,5/115 соот- 
ветствуют ТРДЦН-63000/1 10. На рис. 5 приведены 
графики фазных токов и напряжений в месте воз- 


никновения однофазного короткого замыкания 
(КЗ) К(1). Предшествующий КЗ режим соответ- 
ствует холостому ходу СМ: фазные токи равны О, 
напряжения фаз - номинальные. В момент замы- 
кания напряжение на фазе А упало до нуля и в этой 
же фазе появилась апериодическая составляющая 
и убывающая периодическая. В оставшихся здоро- 
вых фазах ток остался нулевым, а амплитудные 
значения напряжений уменьшились. В токе воз- 
буждения появилась вторая гармоника. Наблюдае- 
мый процесс соответствует классической картине 
для однофазного замыкания [2]. 

В тех же условиях выполнено двухфазное замы- 
кание (рис. 6), фаза А замкнута с фазой В. Волно- 
вые диаграммы соответствуют диаграммам этого 
вида КЗ [2]. 



Рис. 5. Однофазное КЗ на ЛЭП ПО кВ на расстоянии 1 км от шин блочного трансформатора 63 МВА электростанции 
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Рис. 6. Двухфазное КЗ на ЛЭП 110 кВ на расстоянии 1 км от шин блочного трансформатора 63 МВА электростанции 


Проведенные опыты по моделированию элек- 
тромеханических и электромагнитных переходных 
процессов, часть результатов из которых показаны 
на рисунках, выявили полное соответствие поведе- 
ния модели исследуемой системе в классических 
режимах [2, 3, 7]. Это говорит о работоспособности 
цифровой модели, описываемой жесткой системой 


алгебро-дифференциальных уравнений с широким 
диапазоном разброса постоянных времени от 
0,0001 до 30 с, адекватности поведения отдельных 
элементов и всей модели в установившихся и пере- 
ходных режимах реальных электроэнергетических 
систем и ее пригодности к использованию в каче- 
стве учебно-исследовательского инструмента. 
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